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Abstract 

We consider the problem of constraining a particle to a submanifold £ of configuration space 
using a sequence of increasing potentials. We compare the classical and quantum versions of 
this procedure. This leads to new results in both cases: an unbounded energy theorem in the 
classical case, and a quantum averaging theorem. Our two step approach, consisting of an 
expansion in a dilation parameter, followed by averaging in normal directions, emphasizes 
the role of the normal bundle of X, and shows when the limiting phase space will be larger 
(or different) than expected. 



1. Introduction 

Consider a system of non-relativistic particles in a Euclidean configuration space R"+ m whose 
motion is governed by the Hamiltonian 

H=±{p,p) + V(x). (1.1) 
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We are interested in the motion of these particles when their positions are constrained to lie on some 
n-dimensional submanifold £ C K rl+m . In both classical and quantum mechanics there are accepted 
notions about what the constrained motion should be: 

In classical mechanics, the Hamiltonian for the constrained motion is assumed to have the form 
(1.1), but whereas p and x originally denoted variables on the phase space T*R n+m = R"+ m x M"+ m , 
they now are variables on the cotangent bundle T*£. The inner product (p, p) is now computed using 
the metric that £ inherits from M" +m , and V now denotes the restriction of V to £. 

In quantum mechanics, (p,p) is interpreted to mean — A, where A is the Laplace operator, and 
V(x) is the operator of multiplication by V. For unconstrained motion A is the Euclidean Laplacian 
on K" +m , and the Hamiltonian acts in L 2 (R n+m ). For constrained motion, the Laplace operator for £ 
with the inherited metric is used, and the Hilbert space is £ 2 (£, dvol). 

In both cases the description of the constrained motion is intrinsic: it depends only on the Rie- 
mannian structure that £ inherits from K rl+m , but not on other details of the imbedding. 

Of course, a constrained system of particles is an idealization. Instead of particles moving exactly 
on E, one might imagine there is a strong force pushing the particles onto the submanifold. The motion 
of the particles would then be governed by the Hamiltonian 

H x = ±(p,p)+V(x) + \ 4 W(x) (1.2) 

where W is a positive potential vanishing exactly on E and A is large. (The fourth power is just for 
notational convenience later on.) Does the motion described by H\ converge to the intrinsic constrained 
motion as A tends to infinity? Surprisingly, the answer to this question depends on exactly how it is 
asked, and is often no. 

A situation in classical mechanics where the answer is yes is described by Rubin and Ungar 
[RU]. An initial position on E and an initial velocity tangent to E are fixed. Then, for a sequence of A's 
tending to infinity, the subsequent motions under H\ are computed. As A becomes large, these motions 
converge to the intrinsic constrained motion on E. This result is widely known, since it appears in 
Arnold's book [Al] on classical mechanics. However, from the physical point of view, it is neither 
completely natural to require that the initial position lies exactly on E, nor that the initial velocity be 
exactly tangent. Rubin and Ungar also consider what happens if the initial velocity has a component 
in the direction normal to E. In this case, the motion in the normal direction is highly oscillatory, and 
there is an extra potential term, depending on the initial condition, in the Hamiltonian for the limiting 
motion on E. In their proof, E is assumed to have co-dimension one. A more complete result is given 
by Takens [T]. Here the initial conditions are allowed to depend on A in such a way that the initial 
position converges to a point on E and the initial energy remains bounded. (We will give precise 
assumptions below.) Once again, the limiting motion on E is governed by a Hamiltonian with an 
additional potential. Takens noticed that a non-resonance condition on the eigenvalues of the Hessian 
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of the constraining potential W along S is required to prove convergence. He also gave an example 
showing that if the Hessian of W has an eigenvalue crossing, so that the non-resonance condition is 
violated, then there may not be a good notion of limiting motion on E. In his example, he constructs 
two sequences of orbits, each one converging to an orbit on E. These limiting orbits are identical until 
they hit the point on E where the eigenvalues cross. After that, they are different. This means there 
is no differential equation on E governing the limiting motion. For other discussions of the question 
of realizing constraints see [A2] and [G]. A modern survey of the classical mechanical results that 
emphasizes the systematic use of weak convergence is given by Bornemann and Schiitte [BS]. 

The quantum case was considered previously by Tolar [T], da Costa [dCl, dC2] and in the path 
integral literature (see Anderson and Driver [AD]). Related work can also be found in Helffer and 
Sjostrand [HS1] [HS2], who obtained WKB expansions for the ground state, and in Duclos and Exner 
[DE], Figotin and Kuchment [FK], Schatzman [S] and Kuchment and Zeng [KZ]. There are really two 
aspects to the problem of realizing constraints: a large A expansion followed by an averaging procedure 
to deal with highly oscillatory normal motion. Previous work in quantum mechanics concentrated on 
the first aspect (although a related averaging procedure for classical paths with a vanishingly small 
random perturbation can be found in [F] and [FW]). Already a formal large A expansion reveals the 
interesting feature that the limiting Hamiltonian has an extra potential term depending on scalar and 
the mean curvatures. Since the mean curvature is not intrinsic, this potential does depend on the 
imbedding of E in R n+m . 

It is not completely straightforward to formulate a theorem in the quantum case. We have chosen 
a formulation, modeled on the classical mechanical theorems, tracking a sequence of orbits with initial 
positions concentrating on E via dilations in the normal direction. Actually we consider the equivalent 
problem of tracking the evolution of a fixed vector governed by the Hamiltonian H\ conjugated by 
unitary dilations. In order to obtain simple limiting asymptotics for the orbit we must assume that all 
the eigenvalues of the Hessian of the constraining potential W are constant on E. In fact we will assume 
that W is exactly quadratic. Our theorems show that for large A the motion is approximated by the 
motion generated by an averaged limiting Hamiltonian Hb, with superimposed normal oscillations 
generated by X 2 Ho, where Ho is the normal harmonic oscillator Hamiltonian. The Hamiltonians Hb 
and Ho commute, so the motions are independent. These theorems do not require any non-resonance 
conditions on the eigenvalues of the Hessian of W. However, the limiting Hamiltonian Hb does not 
act in L 2 (E), but in L 2 (NT,) where NT is the normal bundle of E. It is only in certain situations where 
one can effectively ignore the motion in the normal directions and obtain a unitary group on L 2 (T) 
implementing the dynamics of the tangential motion. This occurs, for example, if (a) the eigenvalues of 
the Hessian of W are all distinct and non-resonant, (b) the normal bundle is trivial, and (c) we confine 
our attention to a simultaneous eigenspace of all the number operators for the normal motion. In the 
general situation, the dynamics of the additional degrees of freedom in A^E cannot be factored out, 

3 



and we must be content with analysis on L 2 (iVS). 

Our formulation of the quantum theorems invites comparison with the classical mechanical results 
of Rubin and Ungar [RU] and Takens [T]. It turns out that extra potentials that appear in the two cases 
are quite different, and there is no obvious connection. Upon reflection, the reason for this difference is 
clear. If we have a sequence of initial quantum states whose position distribution is being squeezed to 
lie close to S, then by the uncertainty principle, the distribution of initial momenta will be spreading 
out, and thus the initial energy will be unbounded. However, the classical mechanical convergence 
theorems above all deal with bounded energies. The danger in considering unbounded energies is 
that even if the initial energy in the tangential mode is bounded, the coupling between tangential and 
normal modes may result in unbounded tangential energy in finite time. Our assumptions, which allow 
us to obtain a classical theorem despite the unbounded energy, are motivated by quantum mechanics. 
Our results for classical mechanics with unbounded initial energies are quite similar to our results in 
quantum mechanics. 
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2. Classical mechanics: bounded energy 



To give a precise statement of our results we must introduce some notation. The normal bundle 
to £ is the submanifold of M. n+m x W l+m given by 

NT = {(<r,n) : a e T, n £ N a T} 

Here N a T denotes the normal space to S at <j, identified with a subspace of M n+m . 
There is a natural map from TVS into R n + m given by 

t : (a,n) cr + n. 

We now fix a sufficiently small 5 so that this map is a diffeomorphism of NTs = {(cr, n) : \\n\\ < 6} 
onto a tubular neighbourhood of T in R n+m . Then we can pull back the Euclidean metric from M. n+m 
to NTs. Since we are interested in the motion close to T we may use NTs as the classical configuration 
space. This will be convenient in what follows, and is justified below. 

We will want to decompose vectors in the cotangent spaces of NY,$ into horizontal and vertical 
vectors, so we now explain this decomposition. Let tt : NT — > T denote the projection of the normal 
bundle onto the base given by n : (cr, n) cr. The vertical subspace of T ayn NT, is defined to be 
the kernel of dir : T CTi „iVS — ► T a T. The horizontal subspace is then defined to be the orthogonal 
complement (in the pulled back metric) of the vertical subspace. Using the identification of T CTi „7VS 
with T* n NT, given by the metric we obtain a decomposition of cotangent vectors into horizontal and 
vertical components as well. We will denote by (£, rj) the horizontal and vertical components of a 
vector in T^NT. 

The decomposition can be explained more concretely as follows. For each point cr e S, we may 
decompose T cr R" +m = T a T,(BN a T, into the tangent and normal space. Using the natural identification 
of all tangent spaces with M" +m , we may regard this as a decomposition of K" +m . Let Pj and 
be the corresponding orthogonal projections. Since we are thinking of iVE as an n + m-dimensional 
submanifold of R n+m x R n+m , we can identify T( CT , n )NT with the n + m-dimensional subspace of 
R"+ m x W n+m given by all vectors of the form (X, Y) = (<r(0), n(0)), where (a(t),n(t)) is a curve in 
NT, passing through (cr, n) at time t = 0. The inner product of two such tangent vectors is 

((X 1 ,Y 1 ),(X 2 ,Y 2 )) = (X 1 +Y 1 ,X 2 + Y 2 ) (2.1) 

where the inner product on the right is the usual Euclidean inner product. For a tangent vector (X, Y), 
the decomposition into horizontal and vertical vectors is given by 

(X,Y) = (X,PjY) + (0,P?Y) 
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In the statements of our theorems we will want to express the fact that two cotangent vectors, for 
example £\(t) and £(t) in Theorem 2.1, are close, even though they belong to two different cotangent 
spaces. To do this we may use the imbedding to think of the vectors as elements of R 2 (™+ m ). Then it 
makes sense to use the (Euclidean) norm of their difference, \\£\(t) — £(t)\\ to measure how close they 
are. We will use the symbol || • || in this situation, while |£| will denote the norm of £ as a cotangent 
vector. 

We will assume that the constraining potential is a C°° function of the form 

W{a,n) = \{n,A{a)n) (2.2) 

where for each a, A(a) is a positive definite linear transformation on N a T,. The Hamiltonian (1.2) can 
then be written 

H x (a,n,Z,ri) = 1<&0 + ifojj) + V(a + n) + y (n, A(a)n) (2.3) 
Notice that on the boundary of NT,s 1 , for < 5\ < 5, 

H\(a,n,€,r)) > ciA 4 - c 2 

with 

ci = inf W(a + n) > 

(cr,n):<TGS,||n||=5i 

c 2 = sup \V(a + n)\ 

((T,n):ireS, ||n||— <5i 

By conservation of energy, this implies that an orbit under H\ that starts out in NT,s 1 with initial energy 
less than ci A 4 — c 2 can never cross the boundary, and therefore stays in TVE^ . We will only consider 
such orbits in this paper, and therefore are justified in taking our phase space to be T*NT,s, or even 
T*NT, if we extend H\ in some arbitrary way. 

Since we expect the motion in the normal directions to consist of rapid harmonic oscillations, it is 
natural to introduce action variables for this motion. There is one for each distinct eigenvalue ^ (a) 
of A(a). Let P a {a) be the projection onto the eigenspace of u^((t). This projection is defined on N a Y,, 
which we may think of as the range of in R rl+m . Thus the projection is defined on vertical vectors 
in T( CT)n ) NT, and, via the natural identification, on vertical vectors in TT n \NY,. With this notation, the 
corresponding action variable, multiplied by A 2 for notational convenience, is given by 

ry) = ^—( v ,P ari ) + ^±^l( n ,P a n) (2.4) 

Notice that the total normal energy is given by J2 a ^da- The following is a version of the theorem of 
Takens and Rubin, Ungar. 
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Theorem 2.1 Let the Hamiltonian H\ be given by (2.3) where V, W <G C°°, W has the form (2.2) and 
satisfies 

(i) The eigenvalues w 2 (cr) of A(o) have constant multiplicity. 
Suppose that (cr\,n\, £\,r]\) are initial conditions in T*NT,$ satisfying 

(a) \\<r x - cro\\ + Ux - Coll -» 0, 

(b) I*(ax,n x ,Zx,rix)^I°>O t 

as X — > oo. Let ((T\(t),nx(t),t;\(t),r]\(t)) denote the subsequent orbit in T*NT,s under the Hamiltonian H\. 
Suppose that (a(t), is the orbit in T*E with initial conditions (oo, £o) governed by the Hamiltonian 

a 

Then for any T > 

sup \\ax(t)-a(t)\\ + Ux(t)-m\\^0 

0<t<T 

as A — > oo. 

Implicit in this statement is the fact that the approximating orbit stays in the tubular neighbourhood 
for < t < T, provided A is sufficiently large. This theorem is actually true in greater generality. We 
can consider smooth constraining potentials W where \{n, A(a)n) is the first term in an expansion. 
If we choose our tubular neighbourhood so that W(a + n) > c\n\ 2 and impose the non-resonance 
condition u) a (a) ^ ^(a) + w 7 (c) for every choice of a, (3 and 7 and for every a, then the same 
conclusion holds. This theorem is also really a local theorem: if we impose the conditions on W and 
the non-resonance condition locally, and take T to be a number less than the time where a(t) leaves 
the set where condition (i) is true, then the same conclusion holds as well. 

Actually, Takens [T] only treats the case where all the eigenvalues w Q are distinct and the where 
the normal bundle is trivial. On the other hand, he does not require that 7° > 0. This positivity is 
a technical requirement of our proof and arises because action angle co-ordinates are singular on the 
surface = 0. Since Theorem 2.1 is a minor variation of known results, we will not give a proof here. 

3. Classical mechanics: unbounded energy 

We now describe our theorems in classical mechanics where the initial energies are diverging as 
they do in the quantum case. In quantum mechanics, the ground state energy of a harmonic oscillator 
— ^(d/dx) 2 + t;\ 4 uj 2 x 2 is X 2 uj/2. Thus we will assume that the initial values of the action variables J* 
scale like A 2 7°, and therefore that the initial normal energy diverges like A 2 . Examining the effective 
Hamiltonian h (a, £) in Theorem 2.1, one would expect there to be a diverging A 2 >~2 a I®u) a (cr) potential 
term similar to the constraining potential but with strength A 2 . If this potential is not constant, and 
thus has a local minimum (called a mini-well in [HS1, HS2]), no limiting orbit could be expected in 

7 



general unless the initial positions were chosen to converge to such a minimum. For simplicity, we will 
assume that there are no mini-wells, i.e., the frequencies uo a are constant. 

The first step in our analysis is a large A expansion. It is convenient to implement this expansion 
using dilations in the fibre of the normal bundle. It is also convenient to assume that our configuration 
space is all of NY,. This makes no difference, since the orbits we are considering never leave NY,§. 

The dilation d x : NY -> NY is defined by 

d\(a,n) = (a,Xn) 

As with any diffeomorphism of the configuration space, d\ has a symplectic lift D\ to the cotangent 
bundle given by 

D X = d- 1 * = dU 

The expression for D\ in local co-ordinates is given by (5.1). 

Instead of the original Hamiltonian H\ we may now consider the equivalent pulled back Hamil- 
tonian L\ = H\ o D^ 1 . Since D\ is a symplectic transformation, orbits under H\ and orbits under 
L\ are mapped to each other by D\ and its inverse. Therefore, it suffices to study the dynamics of the 
scaled Hamiltonian L\. 

A formal large A expansion yields 

L x = H B + \ 2 H + 0(\- 1 ) 
where Ho is the harmonic oscillator Hamiltonian 

H (<T,n,Z,ri) = ^(v,v) + l(n,A{a)n) (3.1) 
and Hb is the bundle Hamiltonian given by 

H B (a,n,Z,ri) = ±{JZ,Jt) <T + V(v) (3.2) 

The inner product (•, is the inner product on T*£ defined by the imbedding. Here J denotes the 
identification of the horizontal subspace of T* n NY> with the horizontal subspace of T*E given in 
terms of the bundle projection map 7r CTi „ by J — dir*^ . This map is well defined on the horizontal 
subspace, since d7T CTj „ : T a ^ n NY> —> T CT S is an isomorphism when restricted to the horizontal subspace 
of T atn NYi. Thus, its adjoint dir* n is an isomorphism of T*S onto the horizontal subspace of T* n iVE. 
In local co-ordinates Xi , yi defined in section 5 below, where Xi are co-ordinates for S, the map J simply 
identifies dx t e T* jn iVE with dx t E T*E. 

Additional understanding of the Hamiltonians Hb and Hocan be obtained if we introduce another 
metric on NY,. If (X, Y) E T {lJ , n) NY, let 

((X, Y), (X, Y)) x - ||X|| 2 + \- 2 \\P?Y\\ 2 . (3.3) 
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(In Section 7 we describe in what sense this is a limiting form of the pulled-back, scaled, Euclidean 
metric.) If (•, -) A denotes the corresponding metric on the cotangent space, then 

H B + X 2 H Q = V ), (£, r,)) A + y <n, ^(a)n) + V(a) 

The local co-ordinate expressions for Hb and i?o ar e given in (5.9) and (5.10). 

We will use the notation 4>f to denote the Hamiltonian flow governed by the Hamiltonian H. 

Theorem 3.1 Let L\ = H\ o , wfere f/ie Hamiltonian H\ is given by (2.3). Assume that V, W € C°°, 

fe?s the form (2.2), and that the eigenvalues w 2 ofA(a) do not depend on a. 
Suppose that 7a are initial conditions in T*NY, with j\ — > 70 as A — > 00. 
Then for any T > 

sup ^(7 A )-0f- +A2 ^(7 O ) ^0 

0<t<T 

as A — > 00. 

In this theorem the normal energy of the initial conditions, X 2 Ho(j\) grows like A 2 , since Ho(jx) 
is converging to Ho(ja). This leads to increasingly rapid normal oscillations for both orbits <j>f A (7a) 
and (/)f B+A ff °(7o)- Neither orbit converges as A becomes large. It is only their difference that 
converges. 

The convergence of the initial conditions is stated for the scaled variables 7a- To find out what 
this implies for the original variables (<7a, n\, £\, = D7 1 7a we must determine the action of Da on 
horizontal and vertical vectors. This results in the following conditions 

(a) a a — > era, 

(b) An A -> n , 

(c) |a -» >/£o, and 

(d) A _1 ?7a -» »?o 

where (cto, %i £o 5 Vo) = 7o- Here we are thinking of <r, n as vectors in M. n+m and £, 77 as vectors 
in K 2 ("+ m ). We may also compute what these conditions mean for the initial velocities (X\,Y\) e 
Ta\ Six NE, again thought of as vectors in R 2 ("+ m ) . It turns out that 

(c') X x -» X , and 

(d') A-^a - Ko. 

This theorem gives a satisfactory description of the limiting motion if the Poisson bracket of Hb 
and Ho vanishes. Then the flows generated by Hb and Ho commute and the motion is given by the 
rapid oscillations generated by X 2 H Q superimposed on the flow generated by H B - In this situation we 
can perform averaging by simply ignoring the oscillations. 

An example where {Hb, Ho} is zero is when £ has codimension one, or, more generally, if the 
connection form vanishes. Then Hb only involves variables on T*£, so the motion for large A is a 
motion on £ with independent oscillations in the normal variables. The Poisson bracket {Hb, Ho} 



also vanishes if all the frequencies uo a are equal, but in this case the motion generated by Hb need not 
only involve the variables on T*£. 

The motion generated by Hb can be thought of as a generalized minimal coupling type flow. (See 
[GS] for a description of the geometry of this sort of flow.) The flow has the property that the trajectories 
in NY are parallel along their projections onto S. In particular, |n| 2 is preserved by this motion. 

In general, when the frequencies are not all equal, the flows generated by Hb and X 2 Ho interact, 
and Hb + X 2 Ho generates a more complicated flow which need not be simply related to the flows 
generated by Hb and Ho- Let Hb defined by 

H B (l) = lim T- 1 [ H B o <j>?° ( 7 )<ft. (3.4) 

The existence of this limit follows from the Fourier expansion discussed below. This averaged Hamil- 
tonian Poisson commutes with Ho- It turns out that the flow for large A is the one generated by this 
Hamiltonian, with superimposed normal oscillations. 

Theorem 3.2 Assume that the assumptions of Theorem 3.1 hold. Let Ho, Hb and Hb be the Hamiltonians 
given by (3.1), (3.2) and (3.4) respectively. Let 70 G T*NY and T > 0. Then 



sup 

0<t<T 



as A — > 00. 



In this theorem we do not impose a non-resonance condition. However, the form of the averaged 
Hamiltonian H b depends crucially on whether or not resonances are present. 

To explain this further we introduce scaled action variables. Recall that the scaled Hamiltonian 
was defined by L\ = H\ o . We perform a similar scaling on the action variables and define I a by 

I x a oD- x l =\ 2 I a . 



Then 



I a ((T,n,t,ri) = -^-{v,PaV) + ^f{n,P a n). 



Suppose that there are mo distinct eigenvalues lu^ . Then the flows <^[ Q are commuting harmonic 
oscillations in the normal variables. They are periodic, satisfying 4^+2^ = ^° We- therefore obtain a 
group action $ of the m torus T m ° on T*NT, defined by 

for r = (n, . . . , r ma ) e T m °. Notice that 0f = where u = (u u . . . , uj mo ). 
Now we may perform a Fourier expansion of Hb ° & T yielding 

H B o$ T = e i < v ^F v 
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so that 

l/£Z m o 

It turns out that only finitely many F u 's are non-zero. Thus we may exchange the integral and limit in 
the definition of H b with the Fourier sum to obtain 

H B = Y [ lim T_1 / e^^dt) F » = Y F v 
The non-resonance condition on the eigenvalues ui = (wi , . . . , ui mo ) in this situation would be 

Ifz/^Oandi^ ^Othen^w) ^0. (3.6) 

If this condition holds, we find that Hb = Fq. 

We now examine the case mo = m, where there are m distinct frequencies uj a . We wish to describe 
how the limiting motion generated by Hb can be thought of as taking place on E. To begin, since 
{H b, I a } = for each a, each I a is a constant of the motion, so the motion takes place on the level 
sets of Ii, . . . , I m . Furthermore, we want to to disregard the normal oscillations. Technically, we may 
do this by replacing the original phase space T*NY, with its quotient by the group action <f>. This 
amounts to ignoring the angle variables in local action angle co-ordinates. 

It turns out that 

T*7VE/$ = T*£ x R m , (3.7) 

where the variables in E m are the action variables. Since these are constant, we may think of the 
motion as taking place on T*S. To describe the identification (3.7) we first make a new direct sum 
decomposition of each cotangent space T* a n ^NT,. Since there are m distinct eigenvalues uii, . . . , u) m , 
the corresponding eigenvectors, defined globally up to sign, give an orthonormal frame for the normal 
bundle. In this situation the co-ordinates yt = (n, n^a)) are also globally defined up to sign. Thus the 
subspace of n ^ NY, spanned by dy\ , . . . , dy m is globally defined. This subspace is complementary to 
the horizontal subspace, but is not necessarily orthogonal. Given horizontal and vertical components 
(£, rf) of a vector in T* a n ^NY, we may write £ + r\ = £i + r\\ where £i is horizontal and r\\ is in the 
span of dyx, . . . , dy m . The map from T*NY -» T*Y x K m given by 

(o-,n,£,7?) h-> (cr, J£ 1 ,I 1 (a,n,€,r)),...,I m (a,n,£,r))) 

is invariant under <f> and gives rise to the identification (3.7). 

Now suppose that the values of I\ , . . . , I m have been fixed by the initial condition. Then the 
Hamiltonian governing the motion on T*Y depends on these "hidden" variables, and is given by 

h B {u^;h,...,I m ) = \ (£, Oa + V^ + V^h,..., I m ), (3.8) 
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provided the non-resonance condition holds. Given that the eigenvalues are distinct, the following 
implies (3.6) 

lij^k and I ^ m then u)j ± Uk ± u>i ± ui m 7^ (3.9) 

The extra potential V\ is defined in terms of the frame for the normal bundle, n\{a), . . . , n m (er), 
consisting of normalized eigenvectors of A(a). Let bk,i be the associated connection one-form given 
by 

bk,i[-} = (nk,drn[-}) (3.10) 

Then 

V\(a; h,...,I m ) = Y^ I± ^ L \bk,i\ 2 - (3.11) 

k,l UJk 

Notice that the norm \bk.i | is insensitive to the choice of signs for the frame. 

4. Quantum mechanics 

In quantum mechanics, we wish to understand the time evolution generated by H\ for large A, 
where H\ is the Hamiltonian given by (1.2) with (p, p) = — A. As in the classical case, it is convenient 
to replace the original configuration space R"+ m with the normal bundle NT. We will show that if 
the initial conditions in L 2 (M. n+m ) are supported near T then, to a good approximation for large A, the 
time evolution stays near T. Thus we lose nothing by inserting Dirichlet boundary conditions on the 
boundary of the tubular neighbourhood of T, and may transfer our considerations to L 2 (NTg, dvol), 
where dvol is computed using the pulled back metric. If we extend the pulled back metric, and make 
a suitable definition of H\ in the complement of NT,$, we may remove the boundary condition. Thus 
we may assume that that the Hamiltonian H\ acts in L 2 (NT,, dvol). 

More precisely, we let g^s be any complete smooth Riemannian metric on NT, that equals the 
metric induced from the imbedding in the region { (er, n) : ||n|| < e}, for some e < 8. For example, such 
a gNY, could be obtained by smoothly joining the induced metric for small ||n|| with the metric (•, -)i 
given by (3.3) for large ||n||. Let dvol denote the Riemannian density for .gjvs- Let V(<r, n) be a smooth 
bounded function on NT such that V(a 7 n) = V(a + n) when ||n|| < e. Our goal in this section is to 
analyze the time evolution generated by 

H x = ~A + V(a,n) + ^-{n,A(a)n) (4.1) 

acting in L 2 (NT, dvol). Here A denotes the Laplace-Beltrami operator for gNT,- 
We now introduce the group of dilations in the normal directions by defining 



(D x ^)(a,n) = X m l 2 i,(a 1 \n). 
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This is a unitary operator from L 2 (7VS, dvol x ) to L 2 (NT,, dvol) where dvol\ denotes the pulled back 
density dvol\(a, n) = dvo\(<r, A _1 n). Since the spaces L 2 (NT,, dvo\ x ) depend on A, and we want to 
deal with a fixed Hilbert space as A — > oo, we perform an an additional unitary transformation. Let 

(ivoljvs = lim dvo\ x — dvols <8> c^voIr™ 

Then the quotient of densities dvol^ / dvol\ is a function on TVS and we may define M x to be the 
operator of multiplication by ^ 'dvol N ^ / 'dvo\ x . The operator M x is unitary from L 2 (NT,, dvoljvs) to 
L 2 (NT,dvol x ). Let 

U X = D X M X . (4.2) 

Notice that the support of a family of initial conditions of the form U\ip is being squeezed close to £ as 
A — > oo. We want to consider such a sequence of initial conditions. Therefore it is natural to consider 
the conjugated Hamiltonian 

L\ = U* X H X U X , 

since the evolution generated by L x acting on ip is unitarily equivalent to the evolution generated by 
H x acting on U x ip. 

As a first step we perform a large A expansion. Formally, this yields 

L x = H B + \ 2 H o + 0(\- 1 ) 

where Hp is the quantum harmonic oscillator Hamiltonian in the normal variables, and Hb is quantum 
version of the corresponding classical Hamiltonian, except with an additional potential 

K= n{n - 1) S -^\\h\\ 2 . 
4 8 11 11 

Here s is the scalar curvature and h is the mean curvature vector (see equations (7.2) and (7.1)). Notice 
that this extra potential does depend on the imbedding of S in R n+m , since the mean curvature does. 
The quadratic forms for Ho and Hb are 

(4>,H 4>)= [ J(P y #,P v #) ff ,„ + hn,A(a)nM 2 dvol N x (4.3) 

and 

(V,tfsV>=/ \(JP H dip, JP H d^) a + (V(<j,0)+K(a)M 2 dvol N z. (4.4) 

Local co-ordinate expressions for these operators are given by (7.7) and (7.6) below. As in the classical 
case, we can gain additional understanding of these operators by introducing the metric (3.3). Then 

H B + \ 2 H Q = -^A x + y (n, A(a)n) + V(a, 0) + K(a), 

where A x is the Laplace-Beltrami operator on NT, with the metric (3.3). Note that the volume element 
dvolATs is actually \ m times the usual volume element associated to this metric (see Section 7). 
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The operator Hp is explicitly given on C 2 functions in its domain by the formula 




where {nk(cr) : k = 1 . . . to} is any orthonormal basis for NY, and n = Y^ik=i 2/fe n fe(c)- 

It is easy to show that with the metric (3.3), NY is complete so that any positive integer power 
of Hb + A 2 i?o is essentially self-adjoint on for A > [C]. Similarly, because Ho is basically a 
harmonic oscillator Hamiltonian, it is straightforward to show that any positive integer power of Ho 
is essentially self -adjoint on C^°. The operator Hb is more complicated, but also can be shown to be 
essentially self -adjoint on C§°. The argument is not difficult and will be omitted. 

Theorem 4.1 Let gjvs be a complete smooth Riemannian metric on NY that coincides with the induced 
metric zvhen \\n\\ < e,for some e < 6, and suppose V(a, n) is a bounded smooth extension ofV(a + n). Let 
H\ be the Hamiltonian given by (4.1), acting in L 2 (NY, dvol). Assume that A(a) varies smoothly, and that 
the eigenvalues o/cj 2 of A(a) do not depend on a. 

Let L\ = U^H\U\ acting in L 2 (NY, dvobvs). Then, for every ip € L 2 (NY, dvobvs) and every T > 



Just as in the classical case, this theorem provides a satisfactory description of the motion if 
[H B ,Ho] = 0, so that cxp(~it(H B + \ 2 H )) = cxp(-itH B ) exp(-it\ 2 Ho)- As before, this will 
happen, for example, if S has co-dimension one, or if all the frequencies ui a are equal. 

If E has co-dimension one, then the normal bundle is trivial. (We are assuming that £ is compact.) 
Then we have L 2 (NY, rfvoWs) = £ 2 (E,dvol E ) <g> L 2 (R,dy) and H B = h B ® I for a Schrodinger 
operator Hb acting in L 2 (E, dvols). Since Ho = I <B> ho we have that cxp(—it(HB + \ 2 Ho)) = 
exp(-ithB) ®exp{—it\ 2 ho)- This can be interpreted as a motion in L 2 (E, rfvols) with superimposed 
normal oscillations. 

In the case where the frequencies w a are all equal, the normal bundle may be non-trivial, and 
there is not such a simple tensor product decomposition of L 2 (NT,, dvol^s)- However, for some 
initial conditions ip the limiting motion may again be thought of as taking place in i 2 (S, dvols) 
with superimposed oscillations. For example, consider the subspace of functions in L 2 (NY 7 dvoljvs) 
that are radially symmetric in the fibre variable n. This subspace does have a tensor product de- 
composition L 2 (T,,dvo\z) (g> L 2 adial (R m , d m y). It is an invariant subspace for Hb- Furthermore, 
the restriction of Hb to this subspace has the form hs ® I- Thus, if is a radial function in n, 
then cxp(— itL\)ip = exp(— iths) ® cxp(— it\ 2 ho)ipo- As above, we interpret this as motion in 
L 2 (S, (ivols) with superimposed normal oscillations. 

On the other hand, if the normal bundle is non-trivial, it may happen that the limiting motion takes 
place on a space of sections of a vector bundle over S. Instead of giving more details about the general 



lim sup ( e 
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case, we offer the following illustrative example. Instead of a normal bundle, consider the Mobius 
band B defined by R x R / <~, where (x, y) <~ (x + 1, —y). This an 0(1) bundle over S 1 with fibre R. 
An L? function f/ionB can be thought of as a function onlxM satisfying ip(x + 1, — y) = ip(x, y). If 
we decompose ip(x, y), for fixed x, into odd and even functions of y 

^{x, y) = i/> oven 

then tp evcn (x + 1, y) = ip cvcn {x, 2/) and Vodd(z + 1, y) = -ip odd {x, y). (Notice that these are eigen- 
functions for the left regular representation of O(l) on L 2 (R).) Thus Veven can be thought of as an 
L 2 (R, dy) valued function on S 1 , while VWd can be thought of as an L 2 (R, dy) valued section of a line 
bundle over S 1 (which happens to be B itself). In this way we obtain the decomposition 

L 2 (B) = L 2 (S\dx) ® L 2 vcn (R,dy) ® r(S 1 , da) <8> L 2 odd {R, dy) 

where Y is the space of L 2 sections of B. 

In this example, the bundle is flat, so Hb — ~D 2 + V(x) and Ho = ~D 2 acting in L 2 (B, dxdy). 
Leth+ = -D 2 x + V{x) acting in L 2 (S 1 ,dx) and h- = -D 2 x + V{x) acting in T(S 1 , dx). Let h = -D 2 
acting in L 2 (R, dy), with L 2 ven (R, dy) and L 2 dd (R, dy) as invariant subspaces. Then 

e -it(H B + X 2 H ) _ e ~ith + ^ e ~it\ 2 ho ^ e ~ ith - e ~itX 2 h 

So if the initial condition happens to lie in T (g> £ 2 dd , then we would think of the limiting motion as 
taking place in T, with superimposed oscillations in L 2 dd . 

When Hb and Ho do not commute, we perform a quantum version of averaging. Define Hb on 
qf by 

/>oo 

H B i)=\imT~ 1 e itHo H B e- itHo il>dt (4.5) 
It can be shown that H b is essentially self-adjoint. 

Theorem 4.2 Assume that the hypotheses of Theorem 4.1 hold. Let Ho, Hb, and Hb be the Hamiltonians 
defined by (4.3), (4.4) and (4.5). Then, for every ip <G L 2 (NT,, dvolArs) and every T > 

lim sup 

A^oo o<t<T 

The proof that this limit defining H b exists parallels the discussion in classical mechanics. Suppose 
that there are mo distinct eigenvalues uj 2 , . . . , . For each a = 1, . . . , m define the operators I a via 
the quadratic forms 

(V,/ Q V> = / (^-{P V d^,P a P v d^) + ^-(n,P a n)\^\ 2 ) dvobvE 

These operators all commute and satisfy 

a 
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t(H B +\ z H ) 



-it\ z H , 



-itH, 



An expression for I a in terms of local creation and annihilation operators will be given near the end 
of Section 7. In that section we will show that e lTla HBe~ lTl " is periodic in r with period 2tt. Thus if 
we conjugate Hb with e*E r = 7cv , the resulting operator is defined on the torus T m " and has a Fourier 
expansion 

Here r = (n , . . . , r TOo ) and the coefficients F„ are differential operators. As in the classical case, the 
sum is finite. Thus 

e ltH ° H B e~ ltHo = <z it{v,uj) F v . 
This shows that the limit defining H b exists, and is given by 

Hb= F »- 

As in the classical case, we may look for conditions under which the limiting motion can be 
considered to take place on S. Suppose that the eigenvalues u)\, . . . ,uj m are all distinct, and, in 
addition, that the eigenvectors nfe(cr) can be chosen to be smooth functions on all of S. Then the 
normal bundle is trivial, NY = £ x R m and L 2 (NY, rfvoljvs) = L 2 {Y, dvol E ) ® L 2 (R m , d m y). If the 
non-resonance condition (3.9) holds, then 

Hb = (-^ A £ + V^) + K (<r)j ® 1 + 

The term V\ is slightly different from (3.11), because terms arising in its computation do not all commute. 
It is given by 

The joint eigenspaces of I\, . . . , I m are invariant subspaces for H B - The restriction of Hb to such a 
joint eigenspace is the Schrodinger operator — | As + K(a) + if(o-) + Vi, acting in L 2 (S, dvols), where 
Vi is obtained from Vi by replacing the operators Ik by their respective eigenvalues. Thus Hb is a 
direct sum of Schrodinger operators acting in L 2 (S, rfvols). 

5. Co-ordinate expressions 

Our proofs will rely on local co-ordinate expressions for the quantities introduced above. 

Suppose x(a) is a local co-ordinate map for S. Its inverse <j(x) is a local imbedding of M™ onto 
S C M™ +m . Given a local orthonormal frame ni(a), . . . , n m (cr) for the normal bundle, we obtain local 
co-ordinates for NY, by setting 

Xi(<r,n) = Xi(a), i = l,...,n 
yi{cr,ri) = (rii(a),n), i = l,...,m 
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We then may form the standard bases d/dx\, . . . , d/dx n , d/dyi, . . . , d/dy m for the tangent spaces of 
NT, and dx\, . . . , dx n ,dyi, . . . , dy n for the cotangent spaces. This gives rise to local co-ordinates for 
TNT and T*NT in the standard way. For the cotangent bundle, we will denote these by (x, y,p,r) G 
R 2 ("+ m ). Thus (x,y,p,r) denotes the cotangent vector ^pidxi + ^rjdyj in the cotangent space over 

(<KaO>£jyj"i(*))- 

The standard symplectic form for T*NT is the two form given by 

n m 

u) = dpi A dxi + drj A dyj 
i=i j=i 

The dilation map D\ is given in local co-ordinates by 



D x {x,y,p,r) = (x,Xy,p,X V) 



(5.1) 



Clearly this map preserves the symplectic form uj. 

We now compute the local expression for the metric. Let <7i(x) G M. n+m denote the vector 
da(x)/dxi. The tangent vector d/dxi G T^^NT corresponds to the vector in R 2 ("+ m ) given by 
(<7j, J2j Djdnj{(j)[<Ji\). The tangent vector d/dyj corresponds to (0, rij(a)) Here a = a(x), Oi = <Ji(x) 
and n = ^ yjrij(a(x)). Using (2.1) for the inner product, we find that the local expression for the 
metric has block form 

T 



G(x,y) 



B 1 



B~ 




I 


B~ 




'G s + C 0" 




I 


B 


I 







I 




I 







I 



(5.2) 



where Gs = Gs(i) is the metric for S with matrix entries ((Ti(x), aj(x)), B = B(x, y) is the matrix 
with entries 

Bi,j(x,y) = ^2yk{dn k [(Ti\,nj) (5.3) 
k 

and where C = C(x, y) is the matrix with entries 

c i,j( x ,y) = ^2yk({dn k [ui],aj) + ((Ti,dn k [(Jj])) + ^2ykyi(dn k [a i \,dni[a j }) - BB T 

k k,l 

= y^yk((dn k [(Ti],(Tj) + (ai,dn k [(jj]}) + ^2y k yi(dn k [a i ],Pjdni[(T j ]} 
k k,i 

The geometrical meaning of the term Gs + G is given in (7.12) below. 

The inverse can be written 



(5.4) 



G-\x,y) 



'I -B 


T r 


I 





(Gs + G)- 1 
I 



I -B 
/ 



(5.5) 



The local expressions for the projections onto the vertical and horizontal subspaces can now be 
computed. Let Py and Pr denote the projections for the tangent space and P v and P H the projections 
for the cotangent spaces. Then 






0" 




I 


0" 


B T 


J 


Ph = 


-B T 
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and 



GP V G- 



B 
I 



GP H G- 



I -B 




Notice that the vertical subspace of T CTj „7VS is the span of d/dyi, . . . , d/dy m and the horizontal 
subspace of T* n NT, is the span of dx\, . . . , dx n . The map dfn a _ n : T a ^ n NYi — > T CT I] sends d/dxi G 
Tu tn NTi to d/dxi € Ta-Ti and sends d/dyi € T CT „7VE to 0. From this it follows that J = c^r*" 1 , 
defined on the horizontal subspace of T* n NYi sends cfecj € T* n NY, to ctej € T*E. If (<r, n, £, 77) has 
co-ordinates (a;, y,p,r) then £ has co-ordinates 



P"(z,y) 



p 




P - B(x,y)r 


r 








so that J£ has co-ordinates 



p - B{x,y)r. 



We now compute the expressions for ifx, #b and in local co-ordinates. We will abuse notation 
and use the same letters to denote functions on T*iVE and their co-ordinate expressions. Suppose that 
the co-ordinates of (cr, n, £, rj) are (x, y,p,r). Since 



G~ P ^P VT G- i P v = 




J 



we have that 



,G~ L P 



1 &v 



(r, r) 



(5.6) 



(5.7) 



Here, and in what follows, inner products involving co-ordinate vectors always refer to Euclidean 
inner products. For example, (r, r) = Y^T=i r 1- ^ or tne horizontal vectors, we have 



I -B 
I 



P H = P 



H 



so that 



1 uH 



(5.8) 



&0 = (p h \,\ 

= {(p-B^AG^ + CrHp-Br)) 
Therefore the local co-ordinate expression for H\ is 

H x (x,y,p,r) = ^((p- Br) , (Gs + C)^ 1 (p — Br)) + I ( r , r) + y (y, A(x)y) + V(x, y) 



Here C = C(x, y) and £> = B(x, y) are the matrices appearing in the expression for the metric G, 
A(x) is the matrix for A(a) in the basis given by the orthonormal frame n\, . . . n m used to define the 
co-ordinate system and V(x, y) = V(a(x) + 2/fc n fc( cr ( a; ))) ■ Similarly 



H B (x, y,p, r) = -<(p - Br),G^(p - Br)) + V(x, 0) 



(5.9) 
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where B = B(x, y) and Gs = Gx(x). Finally 



H Q {x,y,p,r) = i(r,r) + ^(y,A(x)y) 



(5.10) 



The expressions for Ho and I a simplify if we can choose the vectors in the local orthonormal 
frame to be eigenvectors of A{a). This is always possible if there are no eigenvalue crossings. When, in 
addition, the eigenvalues w 2 (a) do not depend on a there are further simplifications. In what follows 
we will assume that there are mo distinct constant eigenvalues w 2 for a — 1, . . . , m , where uj^ has 
multiplicity k a . We will assume that the local orthonormal frame used to the define the co-ordinate 
system consists of eigenvectors for A(cr). We label them n a j, where a = 1, . . . m and j = 1, . . . , k a 
where for each a, n a j is an eigenvector with eigenvalue L0 2 a . This means that the co-ordinates y and r 
now also acquire a double labelling. 

First of all we have 

H (x, y, p, r) = ^ (r, r) + 1 

a j 

If the co-ordinates of (a, n, £, r?) are (x, y, p, r), then 



(n, P a n) = 



°<,3 



The vertical cotangent vector rj has co-ordinates P l 



. The corresponding tangent vector has co- 



ordinates G~ 1 P V 



which equals 



, by (5.6). Now the projection P a , acting on tangent vectors, 



r 2 ■ 

a, 3 



just picks off the basis vectors d/dy a j, i.e., P a d/dypj = 5p tCc d/dyp t j. Thus 

( V ,P aV )=Y,r 2 

Therefore 

I a (x,y,p,r) = ^J2 r2 c3 + ^-J2'' 



,vl,3 



Notice that in this situation, where the vectors in the local orthonormal frame are eigenvectors of A{a), 
neither Ho nor I a depend on x or p. 

Now we introduce local action-angle co-ordinates. In analogy with creation and destruction 
operators in quantum mechanics, we define the complex quantities 

n . — l£ 



so that 



Va,j — 



' ce,j 



/2uj a 



{a a<j +a* aJ ) 



■ / Ct / 
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The action variables I a j G K and angle variables </? Q j G 5 1 are then defined by 

Notice that ^\ 7 aj j = J a . The change of co-ordinates from (x, y, p, r) to (x, ip, p, I) is symplectic, since 
dr a j A dy a j = dl a j A difaj. This makes it easy to compute the flow §\ a in these co-ordinates. 
Hamilton's equations for the flow are 

±i = 0, pi = 0, i a j = 

<Pa,j = <^3,a 

Thus, under the flow 4>l a each ip a j is translated by t and all the other variables remain unchanged. 
This implies that under the group action $(t), with r = (n, . . . , r mo ) the quantities a QJ evolve as 

G CL qi j ■ 

We now compute the expression for Hb in action angle co-ordinates. We find 

(Br)i = B i,(a,j) {x, v)r a , 

= Yl 6 (aj),(/3,fc)( X ) r a,i2//3,fc 
{3,k,a.j 

V- b \ a .j),((3,k)( X ) , * w * , 

[3 ,k,a.j 

Here bl ~ ,^ ^ (x) — &(a,j),(/3,fe) [ (J i( a; )] is the antisymmetric matrix given by (3.10). The expression for 
Hb is now obtained by substituting this formula for Br into (5.9), which we may rewrite as 

H B (x,p, if, I) = \Y J Vi9 i ' l Vi - J2( Br W' l Pi + \ J2( Br )i9 i,l (Br)i + V(x, 0) 

i,l i,l i,l 

Here g l,i = g % ' l {x) are the matrix elements of (x). To obtain the expression for H B ° 3>(t) we 
simply replace each occurrence of a a .j in the formula above with e ZTa a a j. Since Hb contains only 
constant, quadratic and quartic terms in a a j , a* •, we see that the Fourier expansion of Hb o $(t) has 
finitely many terms, since the v = [u\,..., v ma )'s that appear have J2 a \ v <*\ e {0> 2, 4}. 

6. Proofs of theorems in classical mechanics 

Proof of Theorem 3.1: We begin with some remarks about the co-ordinate charts for T*iV£. We 
will assume that the frames used to defined the co-ordinates consist of eigenvectors of A(a). We 
assume that each chart has the form {(a, n,£,r]) : a £ U,n e N a T,, £ e T* n NH is horizontal, 77 e 
T* n NT, is vertical}, where W is a co-ordinate chart for S. Since S is compact, there is an atlas with 

20 




finitely many charts, and there exists a positive number e\ so that two points in T*NT, both lie in a 
single chart if their projections onto S are a distance less than e\ apart. 
We use the notation 



7a(*) = ^(7a), 7 A (t)^f B+A2ff °(7o). 
Our first estimates are large A bounds on the components of 

7a(*) = fa(t),n A (t), &(*),»?*(*)) 
that follow from the conservation of energy. These bounds are 

\nx{t)\Mx{t)\<C 

and 

!&(*)! < cx 



(6.1) 



(6.2) 



The analogous bounds also hold for 7 x (t) = (a x (t),n x {t),£ x (t),ri x (t)). Clearly |n A (i)| = \y\(t)\ and, 
by (5.7), |r?A(t)| = |?"a(*)|- Thus, (6.1) implies that |j/a(*)| and |?*a (*) j remain bounded. 

To prove these we first consider the action of D x l on^A- Let 7a = (a\, n a, £a, tj\) have co-ordinates 
{x\,y\,p\,r\). Then ^ e T* x nx NT, has co-ordinates 



pH 


P\ 




P\ - B(xx,y\)r x 




r\ 








We now wish to apply D x 1 . Since B(x,y) is linear in y, the scaling in y\ and in r x cancel. In other 
words 

B{x\, X~ 1 y x )Xr x = B(x x ,y x )r x . 

Thus D^ 1 ^ G T* x>x _ lnx NE has the same co-ordinates as £ A G T^^TVE. This implies that as 

A — > 00, 



I^aI 2 = 



Pa - B(x\,y\)r\ 




Pa - B(x\,y\)r\ 




.G-^za.A-^a) 

= ((pa - B(x x ,y x )r x ), (Gy,{x\) + C(xa, A _1 y A )) ^Pa - B(x x ,y x )r x )} 

= IrfTT*- 1 ^! 2 

Thus, for large A, the initial energy satisfies 

£a(7a) = H X o D-Hhx) < llD^txl 2 +C v + y (M 2 + (n x , A(a x )n x )}) 
< CX 2 , 



(6.3) 
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where Cy is an upper bound for V in a neighbourhood of S. Given this bound on the initial energies, 
we may assume that V is bounded, as was explained in the introduction. We now estimate the energy 
for later times t. 

La(7a(*)) = H\ o D-\ lx (t)) > ±\D?b(t)\ 2 Halloo + C\ 2 (|ry A (t)| 2 + \nx(t)\ 2 ) 

> -Halloo + CA 2 (|»7a(*)| 2 + |«a(*)| 2 ) 
Since energy is conserved, i.e., L\("f\(t)) = L\(-f\), this implies (6.1). In a similar way we find that 

\D^Ut)\ 2 <CX 2 . (6.4) 

Now for \y\ < C\ sufficiently large A there is a constant C such that 

G- 1 (x,y)<CG- 1 (x,\- 1 y) 

in any of the finitely many co-ordinate patches. Thus, (6.3) implies 

\Ut)\ < \D^Ut)\, 

so that (6.4) implies (6.2). 

The proof of bounds (6.1) and (6.2) for 7 A (i) is similar. 
We now wish to improve the bound (6.2) to 

\Ut)\ < C (6.5) 

for < t < T. We begin by defining a function Q that depends on our co-ordinate systems. Let 
Xi((j), . . . , Xat(ct) be a partition of unity with each Xk supported in a single co-ordinate patch. Define 
Q = Y1 QkXk, where the local co-ordinate expression for Qk is 

Qk(x,p) - G s (x)-V> + 1. 

(We are abusing notation by using the same letter Qk for the function on T*NT, and its local co-ordinate 
expression.) Given (6.1) we may find a constant C such that 

ICaWI 2 < CQ( 7 a(*)) 

Thus bound (6.5) follows from an upper bound for Q along an orbit. 

To establish such a bound we first estimate the time derivative of Qk (x \ (t) , p\ (t) ) . This derivative 
is given by the Poisson bracket. 

j t Qk(xx(t), Px {t)) ={Q k ,L x }(xx(t), Px (t),px(t),r x (t)) 
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Recall that the orthonormal frame n\ (a), . . . , n m (a) giving our local co-ordinates consists of eigenvec- 
tors of A(a). Thus 

L\ = Hb + \ 2 Hq + E\ 



with 



H B (x,y,p,r) = Q k (x,p) - (B(x, y)r, G^x^p) + \{B(x, y)r, G^^Bix, y)r) + V{x, 0), 

H (x, y, p, r) = ^ (r, r) + ^ ^ J\ yf 

i 

and 

E\{x,y,p,r) = 

\ ((p - B(x, y)r) , ((G E (aO + C(x, X~ 1 y))~ 1 - G^(x)- 1 ) (p - B(x, y)r) ) 
+ V(x,\- 1 y)-V(x,0) 

Since Qk only depends on x and p any Poisson bracket {Qk, F} is given in local co-ordinates by 

8Q k OF 8Q k dF 



{Q k ,F} = J2 



dpi dxi dxi dpi 



Thus {Qk, Ho} = {Qk, Qk} = 0. Using these formulas, together with (6.1) and (6.2) we find 



dt 



Qk{x x {t) lPx {t)) < c (\\ P x(t)\\ 2 + a-IpaWII 



(6.6) 



<CQ k (x x (t),p x (t)) 
Next, writing Hamilton's equations for x\(t) and using (6.1) we find 

dH B 



dp 



\ix(t)\ < 

<CQHx x (t),p x (t)) 

Since the cutoff functions, written in local co-ordinates, only depend on x\ we find that 

\Xk\<C\x x \<CQ~- 



(6.7) 



(6.8 



Now we show if we evaluate Qk and Qj at the same point 7 = (er, n, £ , 77) with \n\ , \r]\ < C then 



\Qk{i)-QAi)\ <CQkh)-i. 



(6.9) 



To see this, we first compute how our co-ordinates change. If (x , y, p, f) are the co-ordinates in the jth 
chart, obtained from the co-ordinates in the ith chart by a change of co-ordinates on £ and a change of 
frame, then 

p = Mp + b 
G^ 1 = M~ 1 G^ 1 M _1 
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where M is the n x n matrix with entries dxi /dxj and 6 is a vector with components r kyid9ki/dxi 
for an orthogonal matrix valued function 9 (x) given by taking inner products of the elements of the 
old and new frames. Thus 



Q fe + 2(6,M- 1 G E 1 p) + l|fc|| 2 + l 



<Q k + CQl 

This implies (6.9). 

Now we are ready to establish a bound for Q along an orbit. Let Q denote dQ(j\(t))/dt. Then 

Q X! QiXi + QjXj 

j 

= QjXj + QiX 3 Xk 

j k,j 

The first term is estimated using (6.6) yielding 

J2Q j Xj<Cj2Q jXj = CQ 

j j 

To estimate the second term, note that since J2k Xk = 1/ we have J2k Xk = 0. Thus 

J^QkXjXk = o 

so that 

XI QjXjXk = J^iQj - Qk)XjXk 

k,j k.j 

< CQ 

by (6.8) and (6.9). Thus we have the differential inequality 

Q<CQ 



which implies 

This implies (6.5) 

Note that (6.7) implies 



Q(lx(t)) < Q(lx(0))e Ct 



\\& x (t)\\,\\a x (t)\\<C (6.10) 



for < t < T. 

We will now show that there exists e > such that if 



lim sup ||7a(t)- 7 A (t)!| =0 (6.11) 

A^oo r g[o,t] 



24 



holds for some t = t\ < T then (6.11) also holds for any t < t\ + e. Since (6.11) holds for t = by the 
assumption on the initial conditions, this will complete the proof. 

So assume that (6.11) holds for t = t\ < T. To compare the two orbits for nearby times, we want 
to ensure that they lie in the same co-ordinate patch. There exists an t\ > such that 7a and 7 A will lie 
in the same co-ordinate chart if \\cr x — a x \\ < e\. 

Choose Ao so that A > Ao implies 

sup || 7 a(t) - 7 V)|| <ei/3 

r£[0,ti] 

Now fix j > j - For t > t\ 

hx{t) - a x (t)\\ < \\a x (t) - a x (h)\\ + \\a x (h) - * x (h)\\ + ||a A (ti) - <7 A (t)|| 
< 2\t-h\C + a/3 

where C is the constant from (6.10). Thus if we choose e < ei/3C then 7a and j x will lie in the same 
co-ordinate chart for t € [ti , t\ + e]. Notice that we do not rule out the the chart changes with A. 

We now write down the differential equation for 7a and 7 A in this common co-ordinate chart. Let 
z e M 2 ("+ m ) denote co-ordinates for T*7V£, i.e., 



z = 



x 

y 
v 

r J 



Denote by z\ the co-ordinates of 7 a and by z x the co-ordinates of 7 A . For a Hamiltonian H, let Xh 
denote the corresponding Hamiltonian vector field given in local co-ordinates by 



X H (z) = 



dH/dx{z) 
dH/dy(z) 
-dH/dp(z) 
-dH/dr(z) 



Then 



d_ 

~di. 



z x {t) = X Hb {z x {t)) + X X 2 Ho (z x (t)) + X Ex (z x (t)) 



(6.12) 



Since Hq is quadratic, the vector field X X 2 Ho is linear, given by 



X X 2 Ho (z) = X 2 Dz 



for a matrix D that is similar to a real antisymmetric matrix. It follows that (6.12) can be written in 
integral form 

z x (t) = e^'-'^ZA^i) + e xHD f e- x ^ D (x Hb (z x (r)) + X Ex (z x (r)))dr 

Jti 
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We may write a similar equation for the co-ordinates of 7 and obtain 

z x (t) - z x (t) = e x2 ^ D (z x (h) - z\t 1 )) 

+e x2tD f e- x2rD (x HB (z x (t))-X Hb (z x (t)) + X Ex (z x (r)))dr 

The harmonic oscillator evolution e x tD is similar to a rotation and therefore uniformly bounded. 
Moreover we have the estimates 

\\X Hb (z x (t)) - X Hb (z x (t)) I) < C \\z x (r) - z x (r)\\ 

and 

\\X Ex (^aM)|| < CA- 1 

These follow from (6.1) and (6.5) which imply that the co-ordinates for the orbits stay in compact sets. 
Thus 

\\z x (t) - z x (t)\\ =C\\z x (t 1 )-z x (h)\\ +C\t-h\ sup \\z x (T)-z x ( T )\\+C\t-t 1 \\- 1 

re[ti,ti+e] 

If we now also insist that 

e < 1/(2C) 

then we find that 

i sup \\z x (r)- z x (t)\\ KCWz^h)- z x (t 1 )\\+Ce\- 1 

1 re[ti,ti+e] 

Since we have only finitely many co-ordinate charts, there is a constant C so that 

C- 1 \\z x (r) - z x (r)\\ < || 7A (r) - 7 A (r)|| < C \\z x (r) z x (r)\\ 
in any chart. Thus we conclude that 

sup || 7 a(t)- 7 A (t)|| ^CW^ihl-^it^W+CeX- 1 

re[ti,ti+e] 



This implies that 



lim sup ||7a(t)- 7 a (t)|| =0 



and completes the proof. 
□ 

Proof of Theorem 3.2: We will show that there exists e > such that if (3.5) holds for some t = t\ <T, 
then (3.5) also holds for any t < t\ + e. So assume that (3.5) holds for some t = t\ <T. 
Define 

^W = ^°°0f B+A2H °(7o) 
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Choosing our co-ordinate charts as in the proof of Theorem 3.1, we find that for small enough e, 
ip\(t) will stay in a single chart for t £ [t\, t\ + e]. This follows from the estimate (6.10) for -f X (t) = 
<f>f B+X H ° (70) and the fact that the harmonic oscillator motion (j)^/ 1 ® keeps the base point a fixed. 

Let w\(t) denote the local co-ordinates of ip\{t). In local co-ordinates, the evolution <f>_ t is 
given by multiplication by e~ tx2 °, and so 

where D is the same matrix, similar to a real antisymmetric matrix, that appeared in the proof of 
Theorem 3.1, and z x (t) are the co-ordinates of r y x (t). Differentiating, we obtain 

dw x (t) _ a - t \ 2 D Y (p t\ 2 D 1 „ r ,u 
— ^ — — e X HB {e w\{t)), 

so that for t g [ii , ti + e], 

/"* 

w x (t)=w x (t 1 ) + e- sX2D X HB (e sX2D w x (s))ds (6.13) 

Now considerthe family ofR 2(n+m ^valuedfunctionson[ti,ti + e]givenby>V = : A > 0}. 

We will show for any sequence Aj — > 00, there is a subsequence Aij suchthatwAu converges uniformly 
to the same limit w M . This will imply that »a — * Wco uniformly. 

The estimates (6.1) and (6.5) of Theorem 3.1 and the fact that the matrices e~ tD are bounded 
uniformly in t imply that W is a bounded family. Moreover, from (6.13) and the boundedness of the 
orbits, it follows that 

\\w x (t)-w x (t')\\<C\t-t'\ 

so that W is equicontinuous. Suppose we are given a sequence Xj — > 00. Then, by Ascoli's theorem, 
there exists subsequence Aij so that w XlJ converges uniformly to w^. We wish to show that is 
always the same, no matter which sequence we start with. Our assumption on t\ implies that w Xl d (t\ ) 
always converges to the same wq, namely to the co-ordinates of <^ B (70). We will show that w^t) is 
the orbit generated by the Hamiltonian H b with initial condition Wo at t = t\. 
Using the uniform boundedness of the matrices e~ tD in (6.13) we find that 

nt 

w oo (t)=w + / e- sX h D X HB {e sX h D Woo {s))ds + o{l) 

Jt! 

as j — > 00. Now e sXl .j D is a symplectic map, being the Hamiltonian flow (f>^° 2 in local co-ordinates. 
It follows that 

e- sX h D X HB (e sX h D Woo (s)) = X H B o<t>»° 2 

1,3 

If we use the Fourier expansion 

H B °<t>" X 2 = z isXl ' i{v,uj) F v 
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we find that 



so that 



H B o<t> " e " 



/t 
e isX l^X F „( Woo {s))ds + o{l) 

Taking j to infinity and using the Riemann-Lebesgue lemma, we find that 



= w Q + X llB (w 00 (s))ds 
Jti 



This identifies w^t) as the orbit generated by H B with initial condition w at t\, as claimed. 
Now we have 



as A — > oo which implies 



sup 

t€[ti,ti+e] 



sup 

te{ti,t!+e] 



e- tx u z\t)~ Woo {t) 







z A (i)-e tA D Woo (t) 



This implies 



sup 

t€[ti,ti+e] 



and completes the proof. □ 



7. More co-ordinate expressions 

In this section we give the co-ordinate expressions that will be needed in our proofs of the quantum 
theorems. 

We begin by defining the second fundamental form, the Weingarten maps and the mean and 
scalar curvatures. Let X and Y be two vector fields tangent to E. Since the Lie bracket [X 7 Y] = 
dY[X] - dX[Y] is tangent to S we find that 

E(X,Y) = P N dX[Y] = P N dY[X] + P N [X,Y] = P N dY[X] 

is symmetric in X and Y. Here P N denotes the projection onto the normal space. By definition, 
II(X,Y) is the second fundamental form. Given an orthonormal frame ni(cr), . . . ,n m (a) for the 
normal bundle, we have 

n(X,Y) = J2(X,S k Y)n k 
k 
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for a collection of symmetric linear transformations Sk on the tangent space. These are called the 
Weingarten maps. Clearly (X, SkY) = (nk,dX[Y]}. But, by differentiating (rife, X) = 0, we obtain 
(dnk[Y],X) + (nk,dX[Y}} = 0, so that the Weingarten maps can also be written as Sk = —P T drik- 
Here P T denotes the orthogonal projection onto the tangent space. 
The mean curvature vector is given by 



h = - Vtr(5 , fe)n fc 

71 ' * 



(7.1) 



fe=i 



while the scalar curvature is 



(7.2) 



Recall that the local expression G(x, y) for the pulled back metric on NT, has the block form (5.2). 
Initially G(x,y) is only defined for ||y|| < S. In our theorem, we wish to extend this metric to a 
complete Riemannian metric on all of NY. One way to achieve this is to join the induced metric for 
small \y\ to the metric (•, -)i given by (3.3) for large \y\. Since the matrix for the metric (•, -)i is 



7 


B 


"G s 


0" 




I 


B 





I 





/ 







I 



the resulting metric on all of NY, would have the matrix 

G(x,y) = 



'I B~ 




'G E +xC 0" 




I B 


/ 




/ 




/ 



i T 



where \ — x{\v\) is a cutoff function that equals 1 for \y\ < e and for |y| > 5. With this special form 
of the extended metric the local co-ordinate expression below remain true on all of NT if C is replaced 
by \C. However, this special form of the extension is not required for our theorems. 



det(Gx 


+ C). Define 






~D X1 ~ 




~ D yi ~ 


D x = 




. D v = 













The local co-ordinate expression for the operator H\ = — \ A + V(a 7 n) + \ 4 W(<r, n) in the region 

|y| < $ is 



D x - BD l 
A, 



n T 



7^2 



(Gs + G)- 1 
/ 



D x - BDy 

A, 



A 4 

+ V(x,y) + —(y,A{x)y) 



\g-V* ((D x - BA) T 5 1/2 (G S + C)~\D X BD v ) + D T y9 ^D v ) 



+ V(x,y) + —(y,A(x)y) 
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Local expressions for the densities on NY, are 



dvol= Jg~&y)\d n x\\d m y\ 



dvo\ x = yJg{x,y/X)\(Px\\d m y\ 



rfvoWE = \/g{^)\d n x\\d m y\ = \/gv(xj\d n x\\d m y\ 
where gT,(x) = det(Gs(a;)). Thus the multiplication operator M\ appearing in (4.2) is multiplication 

ft 



by f x ^ 4 where 



9t.{x) 



We may now compute the local expression for La- Conjugation by D\ results in every multipli- 
cation by a (possibly matrix valued) function F(x, y) being replaced by multiplication by F(x, y/X), 
and every D y being replaced by XD y . Conjugation by M\ simply puts a multiplication by f x l ^ 4 to 
the right of the operator, and a multiplication by f^ 4 to the left. In a co-ordinate system for a domain 
in NT, of the form {(a, n) : a <ElA,n e NY a } let D = ^ x and G\(x, y) be the scaled and extended 



metric taking into account the scaling of D y as well as y. In other words 



G\{x,y) 



I 
XI 



G(x,y/X) 



I 
XI 



(7.3) 



Then 



Lx = - l -fl /i g(x 1 y/X)- 1 ' 2 D T g(x 1 y/Xfl 2 G- x 1 Df- 1,i + V(x,y/X) + ^(y,A(x)y) 
= -\^ 1/ \fl 1/i D T fl'^^G x \fl' 4 Df x 1/i + V(x,y/X) + ^(y,A(x)y) 
Thus in the region where ||y|| < <5A we may use the explicit form of the metric to obtain 

T 



(7.4) 



1/4 -1/2 



D x - BD V 



A, 



1/2 „l/2 

9s fx 



(Gx + Cx)- 1 
X 2 I 



D x - BDy 



-1/4 



A 2 



(7.5) 



+ V(x,y/X) + —(y,A(x)y), 



where C\(x,y) — C(x,y/X). Note that formally putting fx = 1 above, and replacing C\ by 0, we 
obtain for the first line of (7.5) 



1/2 



"At" 


T 


I 


-B 


T 

1/2 


"G^ 1 




I —B 




Dx 


Dy 







I 


5s 


X 2 I 




/ 




Dy 



which is the Laplace-Beltrami operator for the metric which in local co-ordinates is 

T 



/ B 
/ 



G s 

x- 2 i 



I B 
/ 



This is easily seen to be the matrix for the metric (3.3). This explains part of the origin of the Hb + X 2 Ho ■ 
A more complete analysis (to which we now turn) is necessary to understand the origin of the term 
K(a). 
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Before beginning this, note that the local expressions for Hb and Hq are given by 



H B = ^ (As - B(x,y)D y )*G^(D x - B{x,y)D y ) + K{x) + V(x,0) (7.6) 



and 



H a = \D* y D y + l -(y,A{x)y) 



(7.7) 



Here D* and D* denote the formal adjoints with respect to dvolpf^ given by D* = —g^^D^g^/ 2 , 
D* y = -9v 1,2 DTg¥ = -Dl and B* = g^B^g^ = B T . 

We now wish to perform a large A expansion of L x . To state the error estimates precisely, we 
introduce the notation Ek to denote a smooth function of x and y that vanishes to kth order at 
y = 0, evaluated at (x,y/X). Roughly speaking, Ek behaves like (y/X) k for small y/X. The effect of 
differentiating such an error term is given by 



dEk 

dxi 



= E k 



dE k _ / A- 1 ^-! if k > 1 



% [X^Ea if fc = 



In our theorems we will always assume that the eigenvalues w| of A(a) are constant. If we choose 
the orthonormal frame in the definition of our co-ordinates to consist of eigenvectors of A{a) then 
(n, A{a)n) = J2j w | v] ■ We will make this substitution without further comment below. 

Lemma 7.1 In the region where || y\ < SX, the local expression for L\ can be written 
L X = H B + X 2 H a + (D x - BD y )*Ei(D x - BD y ) + E x . 



Proof: In a co-ordinate system for a domain in NT, of the form { (er, n) : a e U, n e NY,^} let D 
and G\(x,y) be given by (7.3). Setting k\ = (1/4) In fx, we may write (7.4) as 

Lx = \{D - dkxYGlHD - dk x ) + V(x,y/X) + ^ $>^ 2 



D x 
A, 



(7.8) 



where dk\ 



d x k\ 
d y k x 



, dk* x = {dk x ) T , and D* = -g^ 1 ' 2 D T g]! 2 . We further expand (7.8) to obtain 



2 2 



Lx = \D*GJ}D + l -dklG^dk x + l -Y, 9v 1,2 (fe /2 (Gx 1 )^ djkx) + V(x, y/X) + y£s , , 

(7.9) 

If || y || < AJthen 



G x 1 (x,y) = 



I -B(x,y) 
/ 



iT r 



(Gv{x) + C{ X ,y/X))-" 
A 2 / 



/ -B{x,y) 
/ 



(7.10) 
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so that in this region we obtain 



Lx = \{D X - BD y )*G^(xy 1 (D x - BD y ) + ^D* y D y 

+ (D x - BD y )*Ei(D x - BD y ) + E 1 + ^J2 i^M + ^ Vi kxf) 

i 

+ V(x,y/\) + ^^y> 

3 

= H B + X 2 H Q + (D x - BD y )*Ei(D x - BD y ) + E 1 



Here we used (d x — Bd y )E k = E k and dk\ 
The lemma will follow if we can show 

A 



Ex 
\~ 1 E 



, so that (d x — Bd y )k x = E^. 



(7.11) 



This requires a more careful expansion of fx- The first step is to uncover the geometrical meaning 
of the term G E (x) + C(x,y) occurring in the expression (5.2) for the metric. Note that 

{dn k [cri],(Tj) = -(Sk<Ti,(Tj) = -(<Ti,Skcrj) = {(Ji,dn k [<jj\) 

and that 

Mk = G^}[{ai, S k <jj)\ 

is the matrix for the Weingarten map S k in the basis <7i , . . . , a n . Let S be the symmetric operator 
defined by (n, II(X, Y)) = (X, SY). Then S — J2k VkSk, an d the matrix for S in the basis <j\, . . . , a n 
is 

M = M(x,y) =J2vkM k {x) 
k 

A short calculation shows 

G E + C = G E (J - Mf 

Given the block form (5.2) of G and (7.12), we obtain 

fx = 9x/gv = det(G(s,y/A)/det(G E (a;)) 

= det(G s (x)(7- A- 1 M(a;,y)) 2 )/dct(G s (a ; )) 
= Act{I-\- 1 M{x 1 y)f. 



Thus 



^ = iln(/ A 1/2 ) = ilndct(/-A- 1 M) 



-trln(7- A _1 M) 

-^X-hT(M) - ^\- 2 tr(M 2 ) + E 3 



^vMS^-jX 2 ^2ykyitr(S k Si) + E 3 



4 
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(7.12) 



This implies that 

1 



and 



2 

1 



Thus 



i 

= \{(tr(S i ))*-tr(S?))-±{i X (S i ))* + E 1 

= — i — « — s~" " 1 

Thus proves (7.11) and completes the proof □ 

We conclude this section by discussing the expression for Hb in local co-ordinates. We may define 
local annihilation and creation operators, using the co-ordinates y a j defined in Section 5, as 



Then we find 



= J2 K,i°«.i + ^ 



We may also write iis in terms of the annihilation and creation operators. We begin with 

{B(x,y)D y )i= hl a ,j,i3.kD Va , o V0,k- 

Notice that the order of D ya . and yp^ is irrelevant here, since b l is antisymmetric in (a, j) and (/3, k). 
Then we can use 

Vp.k = \J~i^( a 0,k + a* f3M ) 

and substitute the resulting expression in (7.6). The resulting formula expresses Hg as a finite sum of 
terms involving product of 0, 2 or 4 annihilation or creation operators. The identities 

e ltHo a a ^e~ ltHo =e- tu "a OJ e itH ° a* aJ e- itH ° =e itu "al J (7.13) 

lead to a finite sum 

e ltHo H B e- itHo = e it{y,u) F v 
that defines the differential operators F v . 
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Lemma 7.2 For ip e C^(NE), e - UH °<p e V(H B ) and 

e ttHo H B e- ltHo ip= F v (f 

where the operators F v are defined by the sum above. 

Proof: It suffices to prove this for <p e Cq° supported in a single co-ordinate patch, since a general 
tp € C§° can be written as a sum of such functions. Introducing our usual local co-ordinates x and y, 
we find that e~ ltH ° is simply a harmonic oscillator time evolution in the y variables. Hence e~ ltH ° tp 
is in Schwartz space. This implies that e~ ztH ° tp e T>(Hb), and that the expansion of Hg into a sum of 
terms involving products of a a j and a* ^ is valid when applied to e~ ltH ° ip. To complete the proof, it 
remains to show that the identities (7.13) hold when applied to a function tp in Schwartz space. This 
follows from 

j t e UHo a a ^e- UHo V = ie UH ° [H Q , a a j}e- itH ° tp 

= iLu a e ltHo [a* aJ a aJ , a aJ ]e- ltHo ip 
= iu ct e itH °[a* aJ ,a a!j }a aJ e- itH °<p 
= -iuj a e UHo a a ^e~ itH °p. 

□ 



8. Proofs of theorems in quantum mechanics 

We begin with two propositions that allow us to transfer our considerations from M. n+m to the 
normal bundle NY,. Let 

d{x,Y,) = inf{||a;-<7|| : a £ £} 
denote the distance to S in R n + m and let 

U s = {x g R n+m : d(x, S) < 5} 

be the tubular neighbourhood of E that is diffeomorphic to ATE 5. The first proposition shows that the 
time evolution in L 2 (M. n+m ) under H\ is approximately the same for large A as the time evolution in 
L 2 (Us) under the same Hamiltonian, except with Dirichlet boundary conditions. 

Proposition 8.1 Suppose that W, V G C°°(R n+m ) with W > and V bounded below. SupposeW(x) = 
if and only if x e E and that W{x) > wq > 0/or large x. 

Suppose X > 1,V e L 2 (R n+m ), \\ip\\ = 1 and \\H x ip\\ < C x \ 2 , where H x = -\& + V + \ i W. Then, 
given e > there exists C2 such that for allt el 

Hi^e-^VII < C 2 \-\ (8.1) 
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Here denotes multiplication by the characteristic function supported on the region indicated in the paren- 
theses. 

Define H s x be the operator in L 2 (U$) given by H\ with Dirichlet boundary conditions on dUs- Then for 
all t e [0,T]andO < e < 6 

\\F {d < €) e- itH ^ - e-^F^W < C 3 \-^ (8.2) 

Here C 2 depends only on C\ and e and C3 depends only onC\,T and e. 

Remark: The power 1 /4 in (8.2) is not optimal. 

Proof: By the assumption on ip and the Schwarz inequality 

Without loss we may assume that V > 0, so that 

1 



2" v " " (8.3) 
(ip,Wip) < Ci\- 2 

It follows that 

C(e)(V,F (d > e) V) < (i/>,F (d > e) Wil;) < C,\- 2 

which proves (8.1), since e~ ltHx ip satisfies the same hypotheses as ip. 
For < ei < a we will need the estimate 

\\F {ei < d < a) ViP\\<C 4 \i, (8.4) 

where C4 depends only on a, e\ and C\. To prove this, choose a function x <= Cg°(M Tl+m ), < \ < 1, 
which is 1 in a neighbourhood of {x : e\ < d(x, S) < a} and vanishes in a neighbourhood of S. Then 

\\F {ei < d < a) V^\\ = ||F (ei < d < a) V(x^)|| < ||V(xV)ll- 
The Schwarz inequality and integration by parts gives 

||V(x^ll<l|A(x^||^||xV'll i 

so that (8.4) follows from 

||A( X ^)|| < C 5 A 2 (8.5) 
and (8.1). To prove (8.5) letp = -iV and calculate, as forms on x C£° 

H 2 X = Jb| 4 + (V + \ 4 W) 2 + Y,Pj( V + Xiw )Pj - \( AV + X * AW ) ( 8 - 6 ) 

3 
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It follows from (8.6) and the fact that C5 is a core for H x that \i> G T^{H\) and 

ll^ 2 X^|| 2 <||i?A(x^)|| 2 + CA 4 , 

or 

\\\ P 2 x n<^c>? + \\Hxn + \\[\p 2 ,xmi 

The last term can be bounded by (8.3), yielding (8.5). 

Let x be a smooth function which satisfies < x < F(d<e/2) an d x = 1 in a neighbourhood of S. 
Because of (8.1) (which holds at t = 0) it is enough to show 

He^xe - *"^ " XV'II < C\- 1/4 

fort e [0,T]. Let 

^A = e^ X e~ 4 * HA V'-X^- 



Integrating the derivative, we obtain 



t,A = i / e ls ^ (tfjfr - X# aK^ 
Jo 

= / e is ^(Vx-p-(i/2)Ax)e- is ^V^, 
Jo 



and thus 



H0 4 ,a!| 2 = / (e- isH Ut,x,(Vx-p-(imAx)e- isH ^)ds. 
Jo 

Let x = 1 on the support of Vx and \ = in a neighbourhood of E. Then from (8.4) 

||<M| 2 < f \\~Xe' lsHi ^x\mX-pe' lsH ^\\+C)ds 
Jo 

< CA5 / ||^ e —^0 t , A ||ds 
Jo 

Now 

(4>t,x,H 5 x <j> t}X ) < 2{ X e- itH ^,Hixe- itH ^) + 2{M,H{*?I>) 

= (e~ itH ^, {H xX 2 + X 2 H X + (Vx) 2 )e- M ^) + (V, {H xX 2 + X 2 #a + (Vx) 2 )^) 

<CX\ 

by the Schwarz inequality. Thus, following the proof of (8.1), 
so that 

HMI 2 < CAiA- 1 

which gives (8.2). □ 
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Since the subset Us C M. n+m is diffeomorphic to NT S C NT, we may think of H{ = -§A + 
V + A 4 Ty as acting in L 2 (NTs, dvol) with Dirichlet boundary conditions on dNT$, where the volume 
form dvol and the Laplace operator A are computed using the pulled back metric, and V and W are 
now the pull backs of the corresponding functions on Us- We may now extend the metric, and the 
potentials V and W, from NTs to all of NT, as explained in Section 4 above. Recall that the extended 
metric is assumed to be complete, that the extended V is bounded and that W = (n, A(a)n) on all of 
NT. We thus obtain an operator H\ acting in L 2 (NT, dvol). Since the extended metric is complete, 
H\ is essentially self-adjoint on C^°. Then it makes sense to talk about e~ ltHx . 

A proposition analogous to Proposition 8.1 holds in this situation, allowing us to approximate 
the evolution under if* with an evolution under H\. For the purposes of this proposition, it does 
not matter how the extensions are made, as long as the conditions on the potentials hold, and the 
state ip that we use for the comparison satisfies H-Ha^II < C\ 2 . Since the statement and proof of this 
proposition are nearly identical to Proposition 8.1 we omit them. 

Having justified the transfer of our considerations to L 2 (NT, dvoljvxi), we now turn to the proof 
of Theorem 4.1. 

Before beginning, we need some quantum energy bounds. 

Lemma 8.2 Let L\ be as in Theorem 4.1 and L ,a = Hb + X 2 Ho- Let L§\ denote either of these operators 
and Rf\ = (A~ 2 Lja + l) \ Let F 2 = -F(| n |/A<e) be a smooth cutoff to the indicated region. When e < 5, 
this cutoff function is supported in the region of NT where the metric is explicitly defined. Let x(cr) be a cutoff 
with support in a single co-ordinate patch. Then, for small enough e and large A, 

II^Kf II + 11X^2 A,^ 2 H + IIA-^A^H < C (8.7) 

If I is a non-negative integer and a, [3 are multi-indices with I + \a\ + |/3| < 2, then 

\\ X F 2 (n) l {\- 1 D x ) a D0R tx \\ < C. (8.8) 

In addition, if I is a positive integer and \a\ + |/3| < 2, then 

\\ X F 2 {n) l {\-'D x TDeR\f\\<C. (8.9) 

Here (n) = + \n\ 2 . 

Proof: Without loss of generality we can assume that V > 1. Set / = %F 2 . Then / <G Cq° with 
< / < 1. Using (7.8) we see that 

L X > \{D dk x yfG^f(D dk x ) + y 5>& ? - 

3 

In the region where / > we can use (7.9) to obtain 
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X 2 I 
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Using \- 2 R^ 2 {L x + X 2 )R X /2 = 1 we obtain 

ll/AX /2 ll < C (8.10) 
A" 1 \\f(D x - BD y - d x k x + Bdyk x )Rl /2 \\ < C (8.11) 

\\{n)R]! 2 \\<C. (8.12) 

On the support of /, d x k x — Bd y k\ is bounded. Thus, using (8.10) and < C\n\ we obtain 
\~ 1 \\fD x R 1 x 2 \\ < C. This proves (8.7) for R\. The proof for i? ,A is similar. 
Define U by L x = \D*G^D + U. Then, using (7.10) we calculate 

L\f 2 L\ = - A {fD*G- x 1 Dy(fD*G- x 1 D)+D*G- x 1 f 2 UD + {Uf) 2 

+ l -D*G- x 1 [D,.f 2 U]+ l -[Uf 2 ,D*]G- x 1 D 

The last two terms above combine to give a multiplication operator given by a function which is easily 
shown to be bounded below by 

-X 2 F|(1 + A 2 M 2 ) 

where \ and F 2 are like \ an d F 2 , with slightly expanded support. It follows that 

^\\fD*G x ~ 1 DR x \\ 2 + X-\fG- x 1/2 \U\^ 2 DR x \\ 2 + \-\fUR x \\ 2 < 1 + A- 4 ||xF 2 (n)Ai? A || 2 

The right side is bounded by (8.7). From \- 2 \\fUR x \\ < C we obtain ||/(n) 2 i? A || < C, which proves 
(8.8) when I = 2. From 

X- 2 \\fG- 1/2 \U\^ 2 DR x \\<C 

we obtain 

\- 1 \\f(n)(D x -BD y )R x \\ <C 

and 

\\f(n)D y R y \\ < C 

which then gives 

\\f{n)\- 1 D x R x \\<C. 

This proves (8.8) when I = 1. Finally we consider the consequences of A _2 ||/£>*G A " 1 -Di?A|| < C. This 
is equivalent to 

A- 2 || J D*G^ 1 ^/i? A || <C 

since the commutator term can be bounded using (8.7). We thus must examine the operator D*G X 1 D 
acting on functions of compact support in K n+m contained in a domain of the form 9a = {(x, y) : 
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\x\ < r, \y\ < eX} When we rescale y — > Xy and D y — > X~ 1 D y , the operator D*G^D goes over 
to an elliptic operator £7 independent of A operating on functions of compact support in a domain 
6 = {(x, y) : \x\ <r,\y\ < e}. The smooth coefficients of the operator E are bounded in 9. It follows 
that if |ck| + \/3\ < 2 

\\D2Dfy\\ <C\\E1>\\ 
for ?/> with support in 6. When we scale back again this implies 

X- 2 \\D a x {XD y ffR x \\<C 

or 

\\(X- 1 D x ) a D f yfR\\\ < C. 

Again, the commutator term which arises from moving / to the left can be bounded using (8.7). This 
takes care of the case I = in (8.8). We have thus proved (8.8) for R x . The proof for Rq jX is similar. 
We now turn to (8.9). We give the proof for R\. The proof for Ro. x is similar. We first show that 

\\f(n) l R{\\<C. (8.13) 

We write / = ff[ where /i has slightly larger support than / and is of the form hi(x)h 2 (\y\/ X). 
Writing /i (n) = g, we have 

g l R l x = gRxg 1 - 1 ^- 1 + g[g l -\ R X ]R 1 ^ 

= gR^g 1 - 1 ^ 1 + gR^Lx^g 1 - 1 }^ 

= gRxg 1 - 1 ^ 1 + gR x {DlJ 1 {n) 1 - 1 + X~ 1 (D X - BD y )* J 2 {n) l ~ 1 + ^{n)'" 1 )^ 
where J\, J 2 and J3 are bounded functions with support contained in supp/i. Thus, from (8.7) 

||^J|<C|| ff , - 1 ^- 1 ||+C||/ 2 <n) , - 1 i2 , A - 1 || 

where f 2 has slightly larger support than f\. Thus (8.13) follows inductively. 

We now let A a . p denote (A _1 D x ) a £)^ and take A = A a . fi with \a\ + \(3\ < 2. Then 

\\g l AR^<\\[A,g l ]R l + 1 \\ + \\Af 2 g l R l + 1 \\ 

where f 2 has slightly larger support than f\. We have 

\\Af2g l R l x +1 \\ < \\Af 2 R x g l R l x \\ + \\Af 2 [g l ,R x ]R l x \\ 

< \\Af 2 R x \\ ■ \\g l R l x \\ + \\Af 2 R x \\ ■ \\ [g l , X- 2 L x ]R l x +1 \\ 

and 

[A 91] = g-u^-ii^D.VD^ 

M+M<i 
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so that 

\\[A, gi }R l x +1 \\ < HW-i^.Xll- 

M + H<1 

where |<7 7;/Jj ;_i| < C(/3(n))' _1 and where fa has slightly larger support than f 2 . Similarly 

[g l , \- 2 L x ] = J^nY^Dy + J 2 {n) l -\\- 1 D x ) + ^(n)'" 1 
where J\, J 2 and J 3 are bounded functions with support contained in supp/i. Thus 

\\[g l ,\- 2 L x ]R l x +1 \\< Y. IIS7W-iA,Xll 

M + M<1 

where |g 7!A1 ,z_i| < C{fz{n)) 1 ^ 1 . Thus again using induction, the result (8.9) follows. □ 
Proof of Theorem 4.1: Since 

it suffices to show. 

lim sup \(ij,e ltLax e- ULx iP} - (tp,tp)\ = (8.14) 

A^oo < t < T 

for a dense set of ip in L 2 (7V£, dvohvs). Let ^ e C§°. Our goal is to show (8.14). 
As a first step, we insert an energy cutoff. Since HLjjaV'II < CA 2 we have 

\\F(L tx /y>^\\ = WF^/x^L^W ■ \\L iXl p\\ 



Set 



- F (Ltx/X 2< ll) 



Then it suffices to show that for each fixed [i > 

lim sup \(F O iTp,e itLo> -e- ltLx F 1 iP}-(F ol TP,F 1 TP)\=0. (8.15) 

A^oo o<t<T 

We now need to show the quantum analogue of the fact in classical mechanics that the orbits stay 
in a bounded region of phase space if we watch the system for a time T < 00 which is independent of 
A. Using energy considerations it follows from Lemma 8.2 that (n) and D y are bounded but only that 
D x cannot grow faster than A. We now seek a A independent bound, showing that up to a fixed time 
T, not too much energy can be transferred from normal to tangential modes. In the quantum setting 
the statement 

\\F 2 D xX e' aL ^F n ^\\ < C, (8.16) 
where F 2 is as in Lemma 8.2, will suffice. 
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We will prove this estimate when L$\ = L\, since the other case when L$\ = Lq\\s similar. Let 
{Xk( a )} b e a partition of unity subordinate to a finite cover of co-ordinate charts. In other words, each 
x\ is supported in a single co-ordinate chart, and J2 k Xk = 1- We- may assume that each Xk is a smooth 
function only of a. Define 

k 

where, in each term, D x and x are defined in terms of the co-ordinates for the chart in which Xk is 
supported. We now want to cut Q off to the region where we have explicit expressions for the metric, 
and then add a constant to regain positivity. So let 

Q = F 2 QF 2 + 1 

Notice that Q and Q commute with F 2 , since in local co-ordinates F 2 is a function of y alone. It is not 
difficult to show that both Q and Q are essentially self -adjoint on C§°(NT,). Define 

q(t) = (F^e^Qe-^F^). 

Then (8.7) follows from 

sup{q(t) :ie [0,T]}<C. 

We will prove a differential inequality as in the classical case. We will need further estimates to 
bound the terms which arise when we compute q(t) and to prove an upper bound for g(0). 

Lemma 8.3 Suppose Fx is a smooth cutoff in the energy \~ 2 L\. Then 



(j\ n ) l {^~ 1 D x ) a Df^jD2xjF 2 FiQ~ 1 



/2 



< c 



ifl + \a\ + \j3\ <2and | 7 | = 1. 

Proof: We use the Helffer-Sjostrand formula (see [D]) 



Fi = J g(z)(R x - z)- x dz!\dz 



where we may take g G C^°(M 2 ) with |g(z)|| lmz\ N < Cm for any N. (We are using the fact that 
F 1 (\- 2 L X )=F 1 (R X ) for Fx e Of (0,2). Let A x - {n) a {\- l D x fD^ X with X e C°°(E), supported 
in the jth co-ordinate patch, xXi = Xi> an d let F 2t x be a smooth function of \n\/X with F 2t iF 2 = F 2 . 
Then 

AxDlx^FxQ- 1 ' 2 = AxF^FxD^XjF^- 1 / 2 + A 1 F 2 s[D2xjF 2 ,F 1 }Q- 1 ^ 2 
Using (8.8), the first term is bounded by a constant times 

WAxF^RxW-WDlxjF.Q-^W <C 
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and it is thus sufficient to show 

\\R^[D2 Xj F 2 ,F 1 }\\<C. 
We compute from the Helffer-Sjostrand formula 

WR^mxjF^FJW < C\\[Dl Xj F 2 ,\- 2 Lx]Rx\\ (8.17) 
For our present purposes we can write 

A 2 

L x = (D x - BD y )*E (D x - BD y ) + -^{D%D y + $> 2 yf) + E 

3 

and we thus obtain 

[D2xjF 2 , \- 2 Lx] = X- 1 D2xj(yF 2 -Dy + Dy \7F 2 ) 

+X- 2 [D2x„ (D x - BD y )*E (D x - BD y )]F 2 + \- 2 E Q 

The first term gives a bounded contribution to (8.17) by Lemma 8.2. The second term can be written 

(A- 1 ^, - BD y )*E \-\D x - BD y ) + D* y E \-\D x - BD y ) 

+ \-\D x - BD y )*E D y + \- 2 E (D x - BD y )^ Xl F 2 

+ \- 1 Dl({d x x j ) T E \-\D x - BD y ) + \-\D x - BD v )*E d x x 3 )F 2 

and again this gives a bounded contribution to (8.17) by Lemma 8.2. □ 
We now return to the proof of Theorem 4.1 and calculate 

q(t) = i(e~ itL ^, Fr[L x , Q}F ie - itL ^). 

Let Fi t i be a function of \~ 2 Lx with slightly larger support than F\, so that F\Fi,\ = Fi. We will 
show that 

Fi,i[iL x ,Q]F ltl <CQ (8.18) 

so that 

q(t) < e Ct q(0). 

First consider any term which arises when the cut-off F 2 = F(\ n yx<e) is differentiated. The 
derivative F' 2 has support in a region of the form {(a,n) : Aei < |n| < Ae 2 } so that F^X/ln]) 1 is 
bounded for any I. Thus = (F 2 {X/\n\) 1 ) X~ l \n\ l so that according to Lemma 8.2, (8.9), such a term 
is bounded (and even decays faster that any inverse power of A. Note that such a term occurring in the 
commutator [Lx, Q] appears alongside D X D^ with \a\ + |/3| < 3 but because we have anFi ; i on the 
left and another on the right, (8.9) even allows \a\ + |/3| < 4 and we still obtain faster than any inverse 
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power of A decay.) Since Q contains the constant 1 such terms are harmless and we will ignore them. 
Thus we are left with showing 

F 1A F 2 [iL x ,Q]F 2 F 1A <CQ. (8.19) 

We write 

h k = D* x G^{x)D x 
when the x refers to the fcth co-ordinate patch. Then 

XkhkXk = ^ (xthk + hkxl) + {d x Xk) T G^}d x Xk 

so that 

[Lx,Q] = £ {^[Lx,xl]h k + \h k [L x ,xl] 

+ [L x ,m k ] + ^xl[L\,h k ] + ^[L x ,h k ]xlj 

where m k = {d x Xk) T G^ 1 d x Xk- We must make use of some cancellation which occurs above so we 
write 

J2\[Lx,xl]h k = ^2\[Lx,XkKhk - hj)x 2 j+J2l^^ h ^ 2 j 

k k.j k,j 

and note that the second term on the right vanishes because J2 k X k — !• Thus we obtain 

l L >» Q] = £ \[Lx,xl]{h k - h 3 ) X 2 3 + \x){h k - h 3 )[L x , X l] 

k,j 

+ [L X ,M] + £ \xl\Lx, h k ] + \[Lx, h k ]xl 

k 

where M = J2k m k- 

In the term [L\, XfcK^fc — hj)x 2 we refer all operators to the jth co-ordinate patch. Thus 

h k - hj = D*G^D X - D* X G^D X 

where ~ refers to the fcth co-ordinate system. We obtain (schematically) D x — M T D X + XE\D y where 
MG^M T = G^ 1 . Hence 

h k - hj = {XEiDy + E )D X + \ 2 E 2 D y D y + XE x D y + E . 

After some calculation we find 

£ \[Lx,Xk](hk - hi) X ) + \x){hk - h 3 )[L XlX 2 k ] 

k,j 

= E XjDKXE^y + E )D xXj + Xl Dl{X 2 E 2 D y D y + \E\D y + E ) 

3 

+ Xj (D* y X 3 E 3 DyDy + X 2 E 2 DyDy + XE^y + E Q ) 

(8.20) 
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where Xj <= C°°(S) with suppxj contained in the jth co-ordinate patch. Noticing the presence of F 2 
in (8.19) and using Lemma 8.3 with a = along with (8.9) of Lemma 8.2, we see that the terms in (8.20) 
give a contribution to the left side of (8.19) which is bounded by CQ. 

We can re-expand M. = M(cr) writing M. = J2k M-x\ an d then we find 

[L A) M] = Xk(D* x Eo + XEiDy + E ) 

k 

which is readily handled by Lemma 8.3 and (8.9) of Lemma 8.2. We now expand the terms involving 
[L\ , ftfe]. After some calculation we obtain 

J2lxl[L-\,h k } + ^[L x ,h k ]xl 

k 

= XkD* x [E X D X + XExDy + \E\ + E ) D xXk 
k 

+ XkD* x ((X 2 E 2 + \Ei)D y Dy + \E\D y + XEi + E ) 

k 

+ J^Xk ((X 2 E 2 + XEi)D y D y + (XEi + E )D y + XE\ + E + X~ 1 E ) 
k 

+ J2(XkD* x E 1 D x + x k E q D x ) 
k 

where Xk € C°°(S) has support in the kth co-ordinate patch with XkXk — Xk- These terms are also 
easily handled with a combination of Lemma 8.2, (8.9) and Lemma 8.3. This completes the proof of 
(8.19) and shows 

q(t) < e ct q(0). 

Finally 

g(0) = (F^QF^) 

has A dependence and must be bounded uniformly in A. But this follows from Lemma 8.3 (with 
I = a = ft = 0) and the fact that \\Q 1/2 tp\\ 2 = (tp, Qtp) < oo, independently of A. 

We now return to (8.15). We introduce a stronger cutoff in the n variable by restricting |n|/A s < 1 
where s e (0, 1). Thus let F 3 = i*(| n |/A«<i) be a smooth cutoff the the indicated region. We note that 

||(1 - F^W < A- S ||(l - FOAVHI • IKn)^!! < CX~° 

by (8.7) of Lemma 8.2. Thus it is sufficient to prove 

lim sup \(FQ^,e ttLa ^F z e- ttLx F 1 %}:) - (F Q . 1 il>,F z F 1 %li)\=0 
A ^°° te[o,T] 

By the fundamental theorem of calculus we obtain 

(F 0tl 1>,e itLo *F 3 e- itL *F 1 1>) - (F ^,F 3 F^) 

= i [ (F ^,e isL ^ ([Lo,x,F 3 ]+ F 3 (L a , x - L x )) F^)ds 
Jo 

(8.21) 
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The term [L .x, F 3 ] contains derivatives of F 3 and thus by Lemma 8.2, (8.9) its contribution to (8.21) 
decays faster than any inverse power of A uniformly for t e [0, T). According to Lemma 7.1, on the 
support of F 3 we have 

Lx - L .x = J2 Xk((D x - BD y )*Ei{D x - BD y ) + E^ Xk . 

k 

Thus, aside from terms involving derivatives of F 3 , which again can be handled by Lemma 8.2, (8.9) 
we need only show that 



lim sup 

A ^°°se[o,T] 



Now 



(Fo 4 e 4sL °.^, (F 2Xk {D x - BD y )*F 3 E 1 (D x - BD y )F 2Xk + ^F^F^^) 







H^iFill <CA-! IKn)^!! <CA- X 
so we need only bound the product 

\\{D X - BD v ) Xk F 2 F^e- sL ^n ■ ll^ill ■ \\(D X - BD v ) Xk F 2 F ie ^ sL ^\\. 



By Lemma 8.2, (8.8) 
and by (8.16) 

Finally 



\BD vXk F 2 F hl \\ < C 



E \\D xXk F 2 F u e isL >^\\ < C. 

s€[0,T] 



ll^i II < c\ s /\ = c\ s -\ 

which proves (8.15) and thus completes the proof of the theorem. □ 

Proof of Theorem 4.2: To prove the theorem, it suffices to show that for any tp g Cq°(NT,), 



lim sup 

A^oc o<t<T 



^ e -it(H B +\ 2 H, 



This can be rewritten as 



lim sup 

A^oo o<t<T 



o) _ e -it\ A Ho e -itH. 



where 



^ x = jt\*Ho e -it{H B +\*H )^ 



We will show that for any ip e L 2 (iVS, dvol^vs) 



sup 

0<t<T 



(V t ,A-e- ltWB V,^) 



= 0, 



(8.22) 



(8.23) 



45 



which will imply (8.22). 

This implication follows from the general fact that if ipt, \ converges to some ipt, oo with 



in the sense that 



SUp ||V*,oo|| < C 
0<t<T 



SUp \{lp ti \ - tp t .oo,f)\ = 
0<t<T 



then, for any continuous function ip t from [0, T] into L 2 (NT,, dvoljvx;), 



SUP \(lp t ,\ - i>t,oo,<Pt)\ = o. 
0<t<T 

To see this, pick an orthonormal basis {(f n }- Then 



SUP \(lpt,\-1pt,oo,¥>t}\ 
0<t<T 



< sup 

0<t<T 



N 



+ sup 

0<t<T 



~y](lpt,X - 1pt,oo,<Pn)(<Pn,<Pt) 
n=l 

N 

<C SUp ^ | {lp t , A - 1pt,oo,<Pn)\ + C SUp ||(1 - P N )<fit\\, 



^2 (lpt,\-1pt,oo,<Pn)(<Pn,<Pt) 
n=N+l 



0<t<T 



where Pm denotes the orthogonal projection onto the span of tpi, . ..,<p at. The first term on the right 
tends to zero as A — > oo, by assumption. Hence 

limsup sup \{i> t ,\-ipt,oo,tPt)\<C sup ||(1 - Pjv)<Pt|| 

A^oo 0<t<T 0<t<T 

But {(j> t : f e [0, T]} is compact and 1 — Pjv goes to zero uniformly on compact sets. Therefore the 
right side tends to zero as N — > oo. 

Thus it suffices to prove (8.23), which we will do in two steps. First, we will show that for every 
sequence Xj — ► oo, there exists a subsequence /Zj and a bounded, weakly continuous L 2 (7V£, dvoljvr) 
valued function ^) t ,oo such that 

sup 1(^,-^,00,^)1^0 (8.24) 

0<t<T 

for every ip e L 2 (NT, 7 dvoljvs). Then, to complete the proof, we will show that ^ tj00 is always the 
same, and equal to e~ ltHB tp. 

To take the first step, we begin with a sequence Xj — > oo. Let {</?n} be an orthonormal basis of 
vectors in Cg°(NY,). Define 

Wn,A(*) = (lpt,X,<Pn) 

Then for fixed n, w„ ; a (i) are a family of functions of i € [0, T], uniformly bounded as A — > oo. Still for 
fixed n, this family is equicontinuous, since the derivative is bounded independently of A. This follows 
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from 



<u\\-\\J2 eltx2{ ^ )F ^\\ 

V 

< M\-^2\\F v ip n \\=C n 

The sum over v is finite. Here we used (7.14), and that ip n is in Cfi°(NY<), and therefore in the domain 
oiF v . 

Using Ascoli's theorem, we may now choose a subsequence Xj 1 of Xj so that wi t \j (t) converges 
to some continuous function W\ t00 (t), uniformly in t for t G [0, T] . Then we may choose a subsequence 
Xj 2 of Xj 1 with W2,x j2 {t) converging uniformly to some continuous function u>2,oo(i)- Continuing in 
this way, and then taking a diagonal subsequence, we end up with a subsequence fij such that 



SUp \W ntH (t) 
0<t<T 



for every n. Notice that 



N 



N 



n=l 



n=l 



<II^H 2 



This implies that Y^= 1 \ w n,oo(t)\ 2 < Ij^ll 2 / so that 

V't.oo = y^Wn,oo(*)yn 
n 

is well defined with || V't.oo 1 1 < ||^||- Clearly, for any n, (V>t, w — ^t,oo> <Pn) — * as j — > oo. This implies 
(8.24) 

Now take the second step of identifying ipt,oo- Let <p <G Cq°(N'E). Then 



,<p) = (1>,<p)+i [\A^,e is ^ Ho H B e- is ^ H °ip)ds 
Jo 

rt 



(8.25) 



Jo 

Since ip € Cg°(NY,) the formula (7.14) implies that 
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Thus the second term of (8.25) tends to zero as j '• ^> oo. On the other hand 

lim f t (^ 00 ,e ls ^ Ho H B e- ls ^ Ho i P }d S = lim V f e is »>^ (^ a>00 , F„<p)ds 

= / (ips,oo,H B (p)ds 
Jo 

by the Riemann Lebesgue lemma. Thus, taking j — > oo in (8.25) we obtain 

(ipt,oo,<p) = (ip,f) +i (ips.oo, H B cp) ds. (8.26) 
Jo 

Now let ip be in the domain of H b- Since is a core for H b, we may use an approximation argument 
to replace ip with e~ lsHs (p and HbP with Hb£~ isHb <P in the equation above. We find, using (8.26), 



d , , _,-.i7„ d 



(is ' (it 



di 



= 

Thus (ip s .oo,e~' lsHB tp) is constant. But when s = 0, equation (8.26) implies (ip s ,oo, e~ tsHB <p>) = (ip, (p). 
Thus(e' sfl ' B ips,oo — ip,0) = for every <p in the domain of Hb- This implies e lsH B ipg^ = tp, or 
V's.oo = e~ lsHB ip, and completes the proof. □ 
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